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Abstract 

We obtain in a closed form the 1/N 2 contribution to the free energy of 
the two Hermitian N x N random matrix model with non symmetric quartic 
potential. From this result, we calculate numerically the Yang-Lee zeros of 
the 2D Ising model on dynamical random graphs with the topology of a torus 
up to n = 16 vertices. They are found to be located on the unit circle 
on the complex fugacity plane. In order to include contributions of even 
higher topologies we calculated analytically the nonperturbative (sum over all 

z (h) 

genus) partition function of the model Z n = J^h^o j^sr ^ or * ne special cases of 
N = 1, 2 and graphs with n < 20 vertices. Once again the Yang-Lee zeros are 
shown numerically to lie on the unit circle on the complex fugacity plane. Our 
results thus generalize previous numerical results on random graphs by going 
beyond the planar approximation and strongly indicate that there might be a 
generalization of the Lee- Yang circle theorem for dynamical random graphs. 

Keywords: Yang-Lee zeros, Lee- Yang theorem, Ising model, random ma- 
trix, random surfaces , 2D gravity . 
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1 Introduction 

Yang and Lee have established long ago []IJ that the statistical theory of phase 
transitions is connected with the distribution of the zeros of the grand partition func- 
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tion on the complex fugacity plane. They have proved that, in the thermodynamic 
limit, those zeros circumscribe closed regions on the fugacity plane where physical 
quantities remanin analytic , thus defining different phases of the system. Although 
complex values for the fugacity are non physical, the physical features of the phase 
transition can be obtained from the study of the distribution of the zeros which 
increase in number and tend to pinch the positive real axis in the termodynamic 
limit . 

In [0 Lee and Yang have proved that the complex zeros of the partition function 
of the ferromagnetic Ising model in a complex magnetic field H lie on the unit 
circle in the complex y-plane (y = e~ 2f3H ). This result, known as the circle theorem, 
makes no assumptions about the details of the lattice like its topology , coordination 
number , etc. Their starting point is the partition function of the model on a given 
static lattice which we can assume to be , for instance , a graph G n with n vertices. 
The corresponding partition function is defined as 

Z(G n ) = £ e (^ G «^' + J ^i*') ) (1) 
{<*} 

where (3 = l/T and Gij = 1 for nearest neighbors and Gij = otherwise. We 
have absorbed the factor f3 in the definition of the magnetic field. As usual the 
configuration of the system is determined by specifying (jj = ±1 on each vertex of 
G n . Henceforth we assume that G n is a closed (no external legs) graph with n four- 
legged vertices. It is important to emphasize that G n does not need to be regular 
and one has Gu = 1 for a link that connects a vertex to itself Defining d as the 
number of spins down we can rewrite Z(G n ) as being proportional to a polynomial 
in the fugacity y, 

Z{G n ) = (cy^)- n E P d y d , (2) 

where c = e _2/3 and each of the coefficients Pd corresponds to a partition function 
of the 2D Ising model on G n without magnetic field (H = 0) with d spins down. 
Therefore they are all positive real and satisfy Pd = P n -d- Based on this and other 
properties of Pd Lee- Yang have proved their theorem || assuring that all zeros of 
Z{G n ) lie on the unit circle : \yk\ = 1, (k = 1, 2, ■ ■ • , n). 

Following Kazakov and Boulatov ^ we can define the Ising model on dynam- 
ical random surfaces (2D Gravity) by treating the lattice G n itself as a degree of 
freedom. Its partition function is obtained from the usual Ising model by summing 
also over all G n with n vertices: 

Z n = £ Z(G n ) = £Z(GW). (3) 

{G„} i 

Here i labels the different graphs with n vertices . 

The above sum is not arbitrary. In our case, where we consider four-legged 
vertices, the weights of the graphs are the corresponding combinatorial factors of a 
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4 field theory. In this case the model of [3,4] is known to have a third order phase 
transition at the critical temperature (3 = log2 from a disordered ( high temperature 
) to an ordered (low temperature ) phase. Clearly Z n has the same form of eq.(Q) 
but with new coefficients Pd = J2%Pd{Gn ) Although the property = P n -d is 
still satisfied , this is not a sufficient condition for the Lee- Yang circle theorem to 
hold . In general the zeros of linear combinations of polynomials have a complicated 
relation to the zeros of the original basic polynomials. However , quite surprisingly 
it has been observed numerically (see || |(|) for dynamical planar graphs (spherical 
topology) wiht 1 < n < 14 vertices that the Yang-Lee zeros are located on the unit 
circle. It is instructive to illustrate this point for planar graphs with a low number 
of vertices (n = 2,3,4) where an analytic analysis is simpler. From eq.(|I|) we easily 
obtain the partition functions given in Fig. 1 . The overall factor (a/ 1 / 2 ) - ™ of a G n 
graph can be canceled by adding a constant to the energy in eq.(Q). Since < c < 1, 





Z(Gj') = c 2 y-' (l+2c 4 y + y 2 ) 




Z(G ( 2 2) ) = c- 2 y- 1 (l+2c 2 y + y 2 ) 



Fig. 1. Planar graphs and partition functions for n = 2 vertices. 



it is easy to verify that the Lee- Yang theorem holds for G 2 and G 2 separately. 
This exemplify, in particular, that the theorem does not depend on specific details 
of the graphs. If we want to treat the graphs as an extra degree of freedom we 
have to sum over G { 2 ] and G ( 2 ] in this case. Taking an arbitrary linear combination 
Z 2 (a, b) = aZ(G { 2 1) )+bZ(G { 2 ) ) one can easily verify that for -2 < b/a < -2(^2-1), 
the zeros of Z 2 (a, b) will not belong to the unit circle. However, in the Kazakov and 
Boulatov's model the combinatorial weights of G 2 and G^ correspond (see ref. J7| 
and the planar Z 2 displayed in [5]) to (a, b) = (2, 16), which brings the zeros to the 
unit circle. 

In the next simplest case, n = 3 , we have four planar graphs, as shown in Fig. 
2. Each graph gives rise to a polynomial of the form 

pW = 1 + y 3 + 3(y + y 2 )a(i) ^ (4) 

where the quantities a® are functions of the temperature which satisfy < < 1, 
(i = 1, 2, 3, 4) . It is easy to demonstrate that each P^ has all its roots on the 
unit circle. Clearly, this is not true for a general linear combination J2t=i ki P3 ■ 
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Z(G ( 3 V C V^ [1+3c 4 (y + y 2 )+ y3] 




Z(0;') = c-y 3/2 [l + 3c 2 (y + y 2 ) + y 3 ] 



Z(G ( 3 3) ) = c-y 3/2 [l+(c 2 +2c 4 )(y + y 2 ) + y 3 ] 




Z(G\ 4> ) = c 3 y" 3/2 [1 + (2c 2 + c 4 )(y + y 2 ) + y 3 ] 



Fig. 2. Planar graphs and partition functions for n = 3 vertices. 



Nevertheless, if hi > 0, the linear combination will have the same basic form of 
up to an overall constant, and we will be back to the unit circle. Obviously, 
the combinatorial factors of the respective graphs (see ref. |7j and the planar Z 3 
displayed in [5]): ki = (32/3, 256/3, 64, 128), belong to this subset. 

For n = 4 vertices there are ten planar graphs. In Fig. 3a and 3b we show 
the diagrams and the corresponding partition functions. Each one corresponds to a 
polynomial, which can be written as: 



P«=l + y * + a f>c 2 + (4 - aHc 4 (y + y A ) + 2 



.(0j 



W 2 i («) 4 
a 2 C + a 3 C 



+ 2 



(4 + a? - 3a^ - 2a£V + (2a^ + - 1 - a£ } )d 



y 2 , 



(5) 



where 



< af < 4, < af, af < 3, 
4 + a f _ 3a « _ 2a f > 0, 2a« + af - 1 - a? > . 



(6) 



Analogous to the n = 3 case, each P4 , under conditions (||), possess only unit 
roots but that is a property that cannot be extended to any linear combination 
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(1) , . 

Z(G 4 ) = c 4 y- 2 [l+2(c 2 +c 4 )(y + y 3 ) 
+ 2y 2 (c 2 +c 4 +c 6 ) + y 4 ] 








, (2). An A 1 

Z(G 4 ') = c 4 y 2 [l+4c 4 (y + y 3 ) 
+ 2y 2 (2c 4 +c 8 ) + y 4 ] 


/ 


n 


\ 


Z(G< 3 >) = ZCc' 1 ') 








Z(cf) = c 4 y- 2 [l+4c 4 (y + y 3 ) 

+ 2y 2 (c 2 +c 6 +c 8 ) + y 4 ] 


Cj 
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Z(G®) = c 4 y" 2 [l+4c 2 (y + y 3 ) 
+ 2y 2 (2c 2 +c 4 ) + y 4 ] 



Fig. 3a. First set of diagrams and partition functions for n = 4 vertices. 



J2}=i kiP^ . However, once again, if ki > (i — 1, • • • , 10) the conditions (|) and 
the polynomial form (|5|) would still hold for the linear combination and we end up 
with unit roots again. Needless to say the combinatorial factors of the graphs 
are all positive and therefore the roots will lie on the unit circle. 

It is remarkable that even for larger n the combinatorial factors of the graphs 
apparently, as our numerical results indicate, are such that the Yang-Lee zeros still 
lie on the unit circle. 

For an arbitrary number of vertices we were not able to make an analytical 
analysis of the location of the Yang-Lee zeros for an arbitrary finite temperature 
but , as one might expect from the results on static lattices, the Yang-Lee zeros 
of Z n should be equally distributed around the unit circle at T = and coalesce 
at y = — 1 as T — > oo . As a self-consistency check we prove that this is indeed 
the case directly from the two Hermitian matrix model in the next section. The 
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Z(G ( f) = c 4 y" 2 [l+4c 2 (y + y 3 ) 
+ 2y 2 (c 2 +2c 4 ) + y 4 ] 


o CS)° 


Z(G ( 4 7 ') = c V 2 [l+2(c 2 +c 4 )(y + y 3 ) 
+ 6c 4 y 2 + y 4 ] 




Z(d*) = c-V 2 [l + 4c 4 (y + y 3 ) 
+ 2y 2 (c 4 +2c 6 ) + y 4 ] 




Z(G 4 ) = c 4 y- 2 [l+(c 2 +3c 4 )(y + y 3 ) 




+ 2y 2 (2c 4 +c 6 ) + y 4 ] 




Z(g" 0> ) = c 4 y" 2 [1 + (c 2 + 3c 4 )(y + y 3 ) 
+ 2y 2 (c 2 +2c 6 ) + y 4 ] 



Fig. 3b. Second set of diagrams and partition functions for n = 4 vertices. 



key ingredient in those cases is the decoupling of the two Hermitian random matrix 
model into the one Hermitian random matrix model. In section 3, using orthogonal 
polynomials, we obtain for arbitrary temperature the 1/N 2 contribution to the free 
energy of the two Hermitian matrix model. In section 4 we present the numerical 
results for the Yang-Lee zeros on graphs of torus topology up to n = 16 vertices. 
In section 5 we calculate, for arbitrary temperature, the non-perturbative partition 
function of the Ising model on random lattices, including all topologies. For this 
end, we have to restrict ourselves to the cases of small matrices N — 1, 2 and graphs 
with n < 20 vertices. Numerical results show that even after summing over all 
topologies we still have the Yang-Lee zeros on the unit circle. 
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2 Free Energy at T = and T — > oo for arbitrary 
topology 

By means of two Hermitian matrices X and F of order NxN, each one associated 
with a value of = =Fl , one can show that the free energy |3|. f| 



1 /jp^ e -Tr[X 2 +Y 2 -2 C Xy+i(e"X4+e- 
#(SS c ? H) = - — log I f -p^ e -Tr{X2+Y2 -2cXY] 



(7) 



is a generating function for partition functions Z n of the 2D Ising model in a constant 
magnetic field H calculated on random graphs G n of n vertices of 4 links each 
according to 

r 

-gc 



E{g,c,H) = -J2 



n=l 



1 



(8) 



Notice that Vfi = T>XT>Y is the usual measure || for Hermitian matrices. The 
quadratic terms of the potential in eq.(^) are responsible for the links (propagators) 

between < ++ >, < > and < H — > while the quartic terms stand for the 

two different sites ( vertices ) (+) and (— ) of the lattices (graphs) each one with 
four links. This corresponds to a 4 interaction field theory in which the fields are 
represented by N x iV matrices. 

Taking ratios of the propagators and comparing with the Boltzmann weights of 
Z(G n ) one identifies the temperature: c = e~ 2/3 . For the specific cases of c = 
(T = 0) and c = 1 (T — > oo), the matrices in eq.(|7|) decouple and the free energy 
E(g,c, H) can be obtained for arbitrary topology in terms of the one Hermitian 
matrix model free energy (E\) as follows. Taking c = (T = 0) in eq.(^) we have 



E{g,0,H) = E 1 



,ge 



H 



+ E 1 ( 



ge 



-H 



(9) 



where 



E ^9) = -^ lo S 



/ -Tr 

JVXe 



ax" 

N 



\ 



JVXe 



- Tr 



X 2 



J 



n=l 



Plugging in eq.(||), we obtain the following expansion in the coupling g, 

E( g ,o,H) = -Y.a n (-gT(i + y n )y- n/2 . 



n=l 



The coefficients a n on their turn have a topological expansion 



An 

h=0 



N 2h 



(10) 



;i2) 
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with h = 0, 1, • • •, corresponding respectively to lattices with spherical topology, 
torus topology, etc. The coefficients offi can be calculated iteratively using orthog- 
onal polynomials for the one matrix model given in (10). From ]7|, [!| we have, for 
instance, the sphere and torus contributions: 

Q (Q) 3"(2n-l)! 

Un ~ n!(n + 2)! [L6) 

a W = ^L»-<W) (14) 
°" 24n { (n!) 2 J ' 1 ' 

In the case c — > 1 (T — > oo), the decoupling of the two matrices in eq.(^) is less 
obvious. With the redefinition g = g(l — c 2 ) 2 we obtain, in the limit c — - > 1 , after 
a change of variables in eq. (0) , 

flfoc-l,!*)^ (f[2/ 1/2 + y- 1/2 ] 

- (_ p )n (1+y) n -n/2 
= - J™ L (1 _ C 2 )2ra ' ( 15 ) 

If we compare eq. (§) with eqs. (|IID and (|i~5D, we obtain the partition functions 

Z n (T = 0) = a n {l + y n )y- n ' 2 (16) 

and 

Z n (T^oo) = a n (l+y) n y- n / 2 . (17) 

Notice that we have dropped the factor [c/ (1 + c 2 ) 2 ]™ since it just amounts to a 
redefinition of the energy by a constant. 

We conclude that for arbitrary number of sites the Yang-Lee zeros of the Ising 
model on dynamical random lattices of arbitrary topology^ are homogeneously dis- 
tributed around the unit circle at T = (yk — e l ^ 2k ~ l \ k = 1,2, • ■ ■ , ra) and 
coalesce at the point yk = — 1 as T — > oo. Both cases coincide with the results of 
the model on a static lattice (see , e.g. , ]10|, [TT]| ) . 



3 Free energy on the Torus at Arbitrary T 

One can reduce the 2iV 2 integrals over Xij, Yij to 2N integrals over their eigen- 
values Xi , yi obtaining || : 



4 Though we have only given a„ and oi , the relations (|l^) and ( |l7|) hold for arbitrary topolo- 
gies. 



s 



where 

/oo roc ^ _____ 
/ ]J dx i dy l u i (c,g,H) _~[(x k - x j ){y k - yj) 
-ooJ-oo. =1 k<j 

with 

Wi = exp ^ - I x i + y i - 2cxiyi + —Xi + — — % I ^ . 

We can further decrease the number of degrees of freedom introducing two 
monic polynomials [|], 

Pj(x) = x j + Y, a kX k 

and 

j'-i 

= y 3 ' + __]hy k , 

k=0 

which are orthogonal with respect to the weight eq.(pO|), 

J J dxdy uj(c,g,H) Pi(x)Qj(y) = Sij hi(c, g, H) . 

Using the relations 

_~[(x k - x{)(y k - y{) = detx 3 ' 1 detyl' 1 = det Py-ifa) det Qj-x{y^ , 

k<l 

one derives: 

fi(9,c,H) 1 h (g,c,H) 

l ) Ln fi ( , c ,H) N m h (0,c,H)> 

where fi = hi/h^i . The last term of eq.(|2S|) can be calculated from eq. 

/oo roo 
/ dy e 
-oo J —oo 

= h {0,c, H) 



r°° POO x 2 +y 2_ 2cxy+ ae x 4 + se_ yi 



_ 1 3gcosh(g) n /J_,, 
iV 2(1 -c 2 ) 2 l iV 2 ' 

The ratios /j and the auxiliary quantities ri, Sj, and tj defined below, 

xPi(x) = Pi+i(x) + riPi-x(x) + SiPi- 3 (x) H 

yQi(y) = Qi+\{v) + qiQ%-i{x) + UQ^iy) h — 



can be determined altogether by solving a set of coupled equations obtained from 
integrals of total derivatives involving eq. fl23|) . Since we are going to look at iV — > oo 
it is more useful to present the continuum version of those equations. Introducing 
the notation: 



x = i/N, e = 1/N , f i+k = Nf{x + ke) 
Qi+k = Nq(x + ke) , r i+k = Nr{x + ke) , 
s i+k = N 2 s(x + ke) , t i+k = N 2 t(x + ke) , 

the continuum equations become ||, §]: 

cq(x) = f(x) {1 + 2ug [r(x + e) + r(x) + r(x — e)]} 



(31) 
(32) 
(33) 



(34) 



x 



cf(x) + — — r(x) = 2ug {s(x + 2e) + s(x + e) + s(x) 

+ r(x) [r(x + e) + r(x) + r(x — e)]} 

2g 



cs[x) 



u 



-f(x)f(x - 2e)f(x - e) 



(35) 
(36) 



and three more equations obtained from ( P4}P^ ) by exchanging (w, r, s, g, t) — > 
(l/u, q, t, r, s), where u = e H . Using the expansion 

fix) = fo(x) + ^/i/a(x) + -^h(x) + ■■■ (37) 

and analogous ones for r(x), q(x), s(x),t(x) and their conjugates, we reproduce, col- 
lecting the (1/N)° terms, the planar equation of ||: 



gx 



c 2 z 3 



z 

V + 3 



-c A + 



Bz 



l-z 



2\2 



9(z) 



(38) 



Here B = 2 [cosh(if) — 1] and z(x) = (6g/c)fo(x) . The above equation furnishes 
fo(x) parametrically. By further collecting the {1/N) terms we get 



1/2 



t 



1/2 



Ql/2 = r l/2 

dto 







d,Sr 



dx 1 dx 
Using such results we derive from the {1/N) 2 terms: 



(39) 



fi{z{x)) 



eg 
2 

+ 2 B' z + 3 c l + c 4 + 16 B z + 28 c z z + 14 B c A z + 75 &z 



C -fz{l- z 2 f {4(1 + zf + 2B{B + 4) [(1 + z 2 f + 32^' 



2 Ji 
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- 5 c 4 z 2 + 112 Bz 3 + 16c 2 z 3 + 8Bc 2 z 3 - 233 c 2 ^ 4 

+ 4 c 4 2 4 + 112 Bz 5 - 308 c 2 z 5 - 154 Bc 2 / + 127 c 2 z 6 
+ 28 c 4 z 6 + 16 Bz 7 + 512 c 2 z 7 + 256 B c 2 z 7 + 217 c 2 

- 98 c 4 z 8 - 268 c 2 - 134 B c 2 z 9 - 271 c 2 z 10 + 154 c 4 z 10 



- 16 cV 1 



8 5 c 2 z 11 + 77 c 2 z 12 - 140 c 4 z 12 + 36 c 2 z 13 



+ 18 5c 2 z 13 + 5cV 4 + 76c 4 z 14 -23c 4 z 16 + 3c 4 z 18 } 



x 



l + z) 



l-c 2 (l-z) 4 ] + 2Bz{l + z 2 ))' 



(40) 



Another contribution of order (1/iV) 2 comes from applying the Euler-Maclaurin 
summation formula : 

N „i 1 1 x 



(*) o + G(iV- 3 ) , (41) 



I E F(z/iV) = £ cfe F(s) + ± [F(0) + F(l)] + -L^ 
on the expression fl26|) which gives 

c, H) =-1 In /io(p, c, if) - jT 1 cfe F(ar) cfe - ^ [F(0) + F(l)] 



12N 2 



Fix) 



+ 0{N' 



where 



F(x) = (1-x) In 



fix, 9 = 0) 



(42) 



(43) 



After using (|38|), we obtain 



g(z{x = 1)) = g 
9{z{x = 0)) = 

1 -x = 1 



(44) 



The integral in (^) can be written, after an integration by parts and using (44) 



as 



I(g,c,H) = £ dx (I - x) In f(x) = I (g,c,H) + 1 ^I 1 (g,c,H) + O(N- 4 ), (45) 



where Iq is the planar contribution: 



1 r 



I (g,c,H) = --l(l-x) 2 \nf(x] 



i 1 r z ( l ) dz' , , ...„ . , 

2 9 2 i z (o) z' L V ;J V ; 
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while I\ stands for the torus contribution, 



h(g,c,H) 



fo(x) 



fi(z') 



(47) 



In the above formulas z(l) is the solution of eq.([3^) with x = 1 that vanishes for 
g = and z(0) = since it is the only solution of eq.(^) for x = that is real in 
the whole range of temperatures. After further manipulations we arrive at 



E(g, c, H) = E®{g, c, H) + ^E^(g, c, H) + 0{N~ A ) 



where 
E^(g,c,B) 



1 z 1 
In- + 



2 g g Jo t 



dt 



g(t) 



i 



2g 2 Jo t 



3 3/2 



(46 



(49) 



Where , following the notation of [4] , the integration limit z = z{x = 1) is of course 
a function of g, c and B. This is the known result from Q for the free energy of the 
Ising on random planar lattices (sphere). The next term which is not known in the 
literature is the torus contribution, 



cosh(_£/~) 1 , 
B«( S ,c,B) = SlT A^ + ^ln 

d R(t) 



+ hn--3gR(z) 
12 g 



+6 / dtg(t) 



d t 



3 

9 -Jo 



dtg 2 (t) 



d R(t) 
d t 



(50) 



Here we have define R{z) = f\(z)/cz. Note that the 1/N contribution vanishes in 



Hq), which is in agreement with the expected topological expansion [12 



4 Yang-Lee Zeros on the Torus 



The free energies and E^ depend on the coupling g explicitly and implicitly 
through z. Now we can expand E^ l > around g = using 



z(x = l,g) 



OO 



k=l 



k\ 



(51) 



where = hk(c, B) can be calculated from consecutive derivatives of g(z) = g with 
respect to g at the point g = 0. Moreover, the integral terms in eq.(|50|) may be 
combined according to 



A(*)= g\t)-2gg{t) 



d R(t) 



d t 



(52) 



and the resulting integral expanded in powers of z, 
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OO I 



d k t Mt) 



t=0 



t=0 



(53) 

is a function of g, c 



In eg. (|53"D , z is given by the series in eq. (|5T|) , and <9 t fc A(t) 
and B. 

Using eqs. (|5T|) and (|53|) together with the expansion of R{z) in powers of g, 
and comparing with eq. (§), we obtain the partition functions on random toroidal 
lattices ZP up to n = 16 vertices . The first 8 results are collected in Table 1. 



zP = KV 1/2 (i + z/) 

Z 2 (1) = § cr V 1 [3 + (4 c 2 + 2 c 4 ) y + 3y 2 ] 

= |c- 3 y- 3 / 2 (1 + 3/) [11+ (-ll + i6c 2 + 17c 4 )?/ + 11 y 2 ] 

Z{ 1] = JL c -4 y -2 [ g37 + ( 1416 c 2 + 1932 c 4) ( y + y 3) + ( 936 c 2 + 2708 c 4 
+ 1304c 6 + 74c 8 )y 2 + 837 y 4 ] 

Zip = % c~ 5 y-^l 2 (1 + y) [579 + (-579 + 1140 c 2 + 1755 c 4 ) + y 3 ) 
+ (579 - 525 c 2 + 805 c 4 + 2245 c 6 + 370 c 8 ) y 2 + 579 y 4 ] 

ZP = fc" 6 y- 3 [21411 + (48492 c 2 + 79974 c 4 ) (y + y 5 ) + (23382 c 2 + 119160 c 4 
+ 142434 c 6 + 36189 c 8 ) (y 2 + y 4 ) + (18630 c 2 + 115452 c 4 + 188488 c 6 
+ 92868 c 8 + 12522 c 10 + 260 c 12 ) y 3 + 21411 y 6 ] 

ZP = f c - 7 y- 7 ' 2 (1 + y) [14571 + (-14571 + 37464 c 2 + 64533 c 4 ) (y + j/ 5 ) 
+ (14571 - 20580 c 2 + 34076 c 4 + 143444 c 6 + 47054 c 8 ) (j/ 2 + y 4 ) 
+ (-14571 + 32319 c 2 + 56581 c 4 + 61278 c 6 + 111853 c 8 + 41307 c 10 
+ 2653 c 12 ) j/ 3 + 14571 j/ 6 ] 

4 1} = § c" 8 jT 4 [710991 + (2050416 c 2 + 3637512 c 4 ) (y + y 7 ) + (884736 c 2 
+ 5729832 c 4 + 9564480 c 6 + 3728700 c 8 ) (y 2 + y 6 ) + (564192 c 2 
+ 5038200 c 4 + 14023200 c 6 + 14328000 c 8 + 5312016 c 10 
+ 549888 c 12 ) (j/ 3 + y 5 ) + (487872 c 2 + 4742568 c 4 + 14910912 c 6 
+ 18385900 c 8 + 9308480 c 10 + 1819608 c 12 + 112832 c 14 + 1198 c 16 ) y 4 
+ 710991 y 8 } 



Table 1: The first eight partition functions on Graphs with torus topology. 
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Fig. 4. Yang-Lee zeros of Z 16 
0.25(+), 0.45(o) and 0.65(x). 



on torus, for temperatures c 



The figure 4 shows the flux of zeros on the unit circle, as the temperature c 
decreases to the critical one, c = 0.25. In Fig. 5 we clearly see that the zeros tend to 
pinch the positive real y— axis as the number of sites increases (thermodynamic limit) 
at the critical temperature c = 0.25. It is worth mentioning that the polynomials 
are different from their planar counterparts, but their zeros still lie on the unit circle. 
The topology slightly changes the position of the zeros along the unit circle. 



5 Non-perturbative Partition Function for Small 
Matrices 

As we have seen in section 3 the exact free energy for spherical and torus contri- 
butions were given as integral representations (eqs. (|49|) and (|5Cf)), in the large N 
limit. These turn out to be the first two contributions of a topological expansion in 
N~ 2h where h stands for the number of handles of the corresponding surface 
in the continuum limit. On the other hand, if we consider small values for N, it 
is possible to obtain an explicit closed form for the free energy (|i~8D including all 
topologies which contribute to a given number of vertices. Perturbatively in g , we 
have in those cases few gaussian integrals to perform ( see (jT9|)). 

It is convenient to write eq. fll~9|) as follows, 



Z(g,c,H) = Z(0,c,H) + Z(g,c,H) 



(54) 
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Fig. 5. Yang-Lee zeros of Z)^ 1 at critical temperature c 
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0.25 for different lattice 



From eqs. (|54] ) and (|T8|) and, after expanding the logarithm, we obtain the free 
energy 



1 



im+l 



m=l 



rn 



Z(g,c,H) 
Z(0,c,H) 



The above expansion can be rewritten in powers of g, 

oo 

Ei p) (c,^) <? 



Z(g,c,H) ^ (p) 



Z(0,c,H) 



(55) 



(56) 



The coefficients be easily calculated for small number of eigenvalues. Here 

we restrict ourselves , for sake of simplicity , to N — 1 (no contribution from the 
vandermonde ) and N = 2 .After rearranging (55) in powers of the coupling g and 
comparing with the expression (8) we obtain the nonperturbative partition functions 
Z n for N — 1,2. We have obtained Z n explicitly for n < 20 sites . In tables 2 and 
3 we display the nonperturbative partition functions Z n (n = 1,2, ...,6) for N — 1 
and 2, respectively. Finally, in Fig. 6 we show the Yang-Lee zeros for our largest 
lattice, n = 20 sites, for N = 2, which are quite close to the ones for N — 1. Here, 
we clearly see the flux of zeros on the unit circle, as the temperature c decreases to 
the critical one, c = 0.25. 
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Zo 


= |c~ 2 w-M2+ (3c 2 + c 4 )v+ 2v 2 l 
2 y l 1 v 1 / y 1 y j 


z. 


= \c- 3 y~ 3 ' 2 (l + y) [11 + (-11 + 16c 2 + 17 c 4 )y + 11 y 2 ] 


z± 


= ! c" 4 y~ 2 [136 + (198 c 2 + 346 c 4 ) (y + y 3 ) + (128 c 2 + 429 c 4 + 
242 c 6 + 17 c 8 )?/ 2 + 136 y 4 } 




— 27 -5 u -5/2 zi , [3714 + f— 3714 + 5440c 2 + 13130 c 4 ) (v + v 3 ) 
+ (3714 - 2800 c 2 + 1455 c 4 + 16370 c 6 + 3545 c 8 ) y 2 + 3714 y 4 ] 


z 6 


= | c" 6 y" 3 [22688 + (33426 c 2 + 102702 c 4 ) (y + y 5 ) + (13056 c 2 + 99855 c 4 
+ 167082 c 6 + 60327 c 8 ) {y 2 + y 4 ) + (9801 c 2 + 89211 c 4 + 199254 c 6 
+ 131184 c 8 + 23601 c 10 + 709 c 12 ) y 3 + 22688 y 6 ] 




Table 2: The first six all topologies partition functions Z n for N = 1. 


6 


Conclusion 



We have obtained a closed and exact expression for the 1/iV 2 contribution to 
the free energy of the two Hermitian random matrix model. From this expression 
one can obtain the partition function of the Ising model in the presence of a 
constant magnetic field H on random graphs (2D Gravity) of torus topology with n 
sites. For arbitrary temperature we have obtained explicitly the partition functions 
Z^ with n < 16. For the specific temperatures c = 0.25, 0.45 and 0.65, we have 
checked that the zeros of the polynomials Z^ (n < 16) lie on the unit circle 
( \y\ = 1) on the complex fugacity plane. 

For the special cases of small matrices lxl and 2 x 2 we have obtained the non- 
perturbative partition function Z n which includes contributions of all topologies. 
We have explicit results for n < 20 sites and again the Yang-Lee zeros lie on the 
unit circle. That is quite surprisingly since there is no Lee- Yang circle theorem for 
dynamical lattices. Taking linear combinations of polynomials in general originates 
non easily predictable changes in their roots. A similar result has been observed 
before for the Z^ [n < 14) on dynamical planar graphs. Our calculations on the 
torus and for higher topologies is thus a strong evidence that the topology of the 
graph plays no special role what the position of the Yang-Lee zeros is concerned. 
The same happens on a static lattice since the Lee- Yang theorem (see 0) is known 
to be independent on the details of the lattice like its topology, number of nearest 
neighbors, etc. It should be stressed that our dynamical lattice consists of sums of 
<p 4 vacuum to vacuum diagrams where the weight of each diagram corresponds to 
its combinatorial factor. Nevertheless, those specific weights seem to be immaterial 
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zi= T^y' 1 ' 2 (i + y) 

Z 2 = ^cr 2 ?/" 1 [17 + (28 c 2 + 6 c 4 )y + 17 y 2 ] 

Z,= |c- 3 y- 3 / 2 (l+y) [7+ (- 7 + 12c 2 + 9c 4 )y + 7y 2 ] 

Z 4 = c" 4 y- 2 [2011 + (3648 c 2 + 4396 c 4 ) (y + y 3 ) + (2528 c 2 + 6424 c 4 
+ 2952 c 6 + 162 c 8 ) y 2 + 2011 y 4 ] 

Z 5 = y|oC" 5 ?/" 5/2 (1 + y) [11394+ (-11394 + 21740 c 2 + 35230 c 4 ) (y + y 3 ) 
+ (11394 - 9800 c 2 + 14055 c 4 + 44770 c 6 + 7945 c 8 ) y 2 + 11394y 4 ] 

Z 6 = c~ 6 y~ 3 [252169 + (501588 c 2 + 1011426 c 4 ) (y + y 5 ) + (231528 c 2 

+ 1299915 c 4 + 1737966 c 6 + 513126 c 8 ) (y 2 + y 4 ) + (183438 c 2 + 1229793 c 4 
+2214752 c 6 + 1222392 c 8 + 188238 c 10 + 4767 c 12 ) y 3 + 252169 y 6 ] 



Table 3: The first six all topologies partition functions for N = 2. 

concerning the unit circle. It is remarkable that replacing those specific weights 
by arbitrary positive constants, at least for n < 4 sites, we still have the roots 
on the unit circle. That happens on the sphere, on the torus and also for linear 
combinations, with positive constants, of both topologies. In fact this is probably 
the reason for the location on the unit circle of the zeros of the nonperturbative 
partition functions that we have obtained for the special cases of 1 x 1 and 2x2 
matrices. In those cases we have a linear combination of all topologies wiht the 
positive coefficiens N~ 2h where h — 0, 1, • • • oo labels the different topologies . Such 
results lead us to conjecture that linear combinations , with positive coefficients , of 
partition functions of the 2D Ising model calculated on different graphs G n have all 
their n Yang-Lee zeros on the unit circle. At this point a word of caution is in order 
, namely , Lee and Yang have proven their circle theorem for a class of polynomials 
which includes the partition function of the Ising Model in a constant magnetic field 
, nevertheless , it is easy to find numerical examples of linear combinations of such 
polynomials which donot obey the circle theorem even if we take positive coefficients. 
Therefore it is not obvious how the Lee- Yang theorem can be generalized for sums 
over the graphs (lattices). It cannot be discarded that the polynomials we have for 
the dynamical case correspond to a new class of polynomials not encompassed by the 
Lee- Yang theorem but whose roots still lie on the unit circle. A clear understanding 
of this issue demands more work which is now in progress. We emphasize that we 
have only looked at the global location of the Yang-Lee zeros rather than their local 
distribuition close to the positive real axis and its relationship with the topology. 
This is now under investigation. Finally, it is worth pointing out that we were able 
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Fig. 6. Yang-Lee zeros of Z20 with N = 2, for temperatures c = 
0.25(+), 0.45(o) and 0.65(x). 



to sum over all topologies for a given finite number of vertices because the sum is 
actually finite since not all topologies can contribute to a graph with a finite number 
of vertices with finite number of links. Thus , we have a natural cut-off . 



7 Note Added 



After finishing this work we became aware of |14| where for the special case N = 1 



("Thin Graphs") it was shown in the thermodynamic limit that the closest zero to 
the positive real axis is located on the unit circle on the complex fugacity plane. 
That is clearly in agreement with the expectations from our finite size calculations 
of all Yang-Lee zeros in the case N — 1 (see section 5). 



8 Acknowledgements 

D.D. is indebted to R.Schrock for an interesting discussion. This work was partially 
supported by FAPESP under contracts # 96/11747-2 , # 97/06499-2 , # 98/12626- 
(D.D.) and # 97/01378-2 (L.C.A.). L.C.A. would like to thank the Departamento 
de Fisica Matematica of Universidade de Sao Paulo for their hospitality. 



18 



References 

[I] C.N. Yang and T.D. Lee, Phys. Rev. 87 (1952) 404. 
[2] T.D. Lee and C.N. Yang, Phys. Rev. 87 (1952) 410. 
[3] V.A. Kazakov, Phys. Lett. B119 (1986) 140. 

[4] D. Boulatov and V.A. Kazakov Phys. Lett. B186 (1987) 379. 
[5] M. Staudacher, Nucl. Phys. B336 (1990) 349. 

[6] J. Ambjorn, K.N. Anagnostopoulos and U. Magnea, Nucl. Phys. (Proc. Suppl.) 
B63 (1998) 751 ( |hep-lat/9708"0"T4l) , Ibidem Mod. Phys. Lett. A12 (1997) 1605 
( |hep-lat/9705004j ). 

[7] E. Brezin, C. Itzykson, G. Parisi and J.B. Zuber, Commun. Math. Phys. 59 
(1978) 35. 

[8] M.L. Mehta, Commun. Math. Phys. 79 (1981) 327. 

[9] D. Bessis, Commun. Math. Phys. 69 (1979) 147. 

[10] S.-Y. Kim and R.J. Creswick, Phys. Rev. Lett. 81 (1998) 2000. 

[II] V. Matveev and R. Shrock, Phys. Rev. E53 (1996) 254; Ibid., Phys. Lett. 
A215 (1996) 271. 

[12] P. Di Francesco, P. Ginsparg and J. Zinn- Justin, Phys. Reports 254 (1995) 1. 
[13] C. Itzykson and J.B. Zuber, J. Math. Phys. 21 (1980) 411. 
[14] D.A. Johnston , Journal of Phys. A31 (1998) 5641. 



Figure Captions: 

Fig. 1. Diagrams and partition functions for n = 2 vertices. 

Fig. 2. Diagrams and partition functions for n = 3 vertices. 

Fig. 3a. First set of diagrams and partition functions for n = 4 vertices. 

Fig. 3b. Second set of diagrams and partition functions for n = 4 vertices. 

Fig. 4. Yang-Lee zeros of Z$ on torus, for temperatures c = 0.25(+), 0.45(o) and 0.65(x). 

Fig. 5. Yang-Lee zeros of at critical temperature c = 0.25 for different lattice 

sizes n on torus: n = 10(x), 13(o) and 16(+). 

Fig. 6. Yang-Lee zeros of Z 2 o with N = 2, for temperatures c = 0.25(+), 0.45(o) and 0.65(x). 
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